Theory of plastic vortex creep 
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We develop a theory for plastic flux creep in a topologically disordered vortex solid phase in type- 
II superconductors. We propose a detailed description of the plastic vortex creep of the dislocated, 
amorphous vortex glass in terms of motion of dislocations driven by a transport current j. The 
plastic barriers Upi{j) oc show power-law divergence at small drives with exponents n = 1 for 
single dislocation creep and n = 2/5 for creep of dislocation bundles. The suppression of the creep 
rate is a hallmark of the transition from the topologically ordered vortex lattice to an amorphous 
vortex glass, reflecting a jump in fj, from /i — 2/11, characterizing creep in the topologically ordered 
vortex lattice near the transition, to its plastic values. The lower creep rates explain the observed 
increase in apparent critical currents in the dislocated vortex glass. 

PACS numbers: 74.60. Ge, 74.60. Jg, 61.72.Bb 

One of the most fascinating dynamic phenomena of 
complex systems with internal degrees of freedom is the 
thermally activated motion of elastic media in a random 
environment (creep) characterized by a highly nonlinear 
response to a dc driving force, F: v ^ exp (— const /Ti^^), 
where v is the velocity, T is the temperature, and fj, is the 
exponent depending on the geometry and the dimension- 
ality of the driven medium. The concept of thermally ac- 
tivated creep, which was originally introduced to describe 
the driven motion of elastic manifolds through quenched 
disorder in the context of the dynamics of dislocation 
and/or domain walls in inhomogeneous environments 
was extended later to driven vortex lattices and 
has eventually become a paradigm for the thermally ac- 
tivated dynamics of disordered media. The concept of 
creep proved to be especially successful in describing a 
wealth of the low temperature transport and relaxational 
properties of the vortex state of high temperature super- 
conductors 1^. Recent observations j|] on 2D magnets 
confirmed that domain walls also exhibit creep behavior. 
Creep dynamics results from the scaling of energy barri- 
ers U{j) ^ which control thermally activated motion. 
The derivation of such energy barriers diverging at small 
currents (driving forces) was based on the elastic behavior 
of the pinned structures; thus in the common view creep 
behavior is implicitly attributed to the elastic medium 
free of topological defects. 

The description of thermally activated dynamics of 
amorphous structures containing a large amount of topo- 
logical defects remains a major challenge. This is a long 
standing problem in the theory of work hardening and 
related relaxation processes in dislocated solids. In the 
context of vortex physics the quest for the description 
of thermally activated dynamics of a topologically disor- 
dered medium was motivated by the observation of the 
disorder-induced transition between a low-field quasilat- 
tice or Bragg glass (BrG) [| 10 1, the phase free of topo- 
logical defects, and a high-field phase characterized by an 
enhanced apparent critical current The latter phase 
was suggested to be a topologically disordered vortex state. 
It was shown that the low-field BrG phase is indeed stable 



against dislocation formation ||Tl|,|T2|. At high magnetic 
fields the vortex solid undergoes a structural transition 
which is described ||T^-[T5|] as a topological transition be- 
tween the BrG and a dislocated, amorphous VG (AVG) 
jl^Jl^. In a recent series of experiments |^,^ the phase 
coexistence characteristic for a first-order transition was 
established, and creep barriers in the high-field vortex 
state were shown to diverge faster than creep barriers in 
the low-field elastic phase. 

These recent experimental findings call urgently for a 
theoretical description of the AVG, and the main the- 
oretical problem is to find an appropriate quantity en- 
abling parameterization of the amorphous phase. A first 
step towards such a quantitative description was made in 
Ref. where all phase transitions between vortex lattice 
phases were described in terms of dislocation-mediated be- 
havior, and a free energy functional .F[n£)] for an ensemble 
of directed dislocations in the presence of thermal fluctu- 
ations and quenched disorder was derived {n jj is the areal 
dislocation density). The BrG- AVG transition was found 
to be of weak first order in accordance with the exper- 
imental results of Ref. at the BrG- AVG transition 
dislocations enter with a density no — given by the 
positional correlation length Ra on which typical vortex 
displacements are of the order of the lattice spacing a 
§. Upon increasing the magnetic field up to the criti- 
cal point the dislocation density of the AVG increases to 



vortex-liquid-like values njj 



such that the AVG 



and vortex liquid phases become thermodynamically in- 
distinguishable at the critical point. 

In this Letter, building on the aforementioned ideas, 
we propose a quantitative description of plastic creep in 
terms of the dislocation degrees of freedom. We establish 
and find a critical plastic current jpi below which dis- 
locations are collectively pinned and plastic creep occurs 
via the activated motion of collectively pinned dislocation 
lines. The critical plastic current is lower than the critical 
current for vortex depinning jpi < jc, hence plastic mo- 
tion of depinned dislocations sets in before viscous flow of 
the entire vortex lattice can occur. We calculate the asso- 
ciated plastic plastic creep energy barriers Upi{j) j^'^''' 
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diverging infinitely at j 0. We treat dislocations in 
the pinned vortex lattice as directed elastic strings sub- 
ject to a pinning force which we derive from the Peach- 
Kohler force exerted on vortices by the pinning centers. 
We show that an external current sent through a dislo- 
cated vortex lattice generates a Peach-Kohler force with 
a driving component for dislocation glide. Knowledge of 
the pinning and driving forces on the elastic dislocation 
string enables us finally to study the glassy dynamics of 
the dislocation, in particular the depinning threshold for 
dislocation glide and the energy barriers for plastic creep 
below the depinning threshold. 

The energy of a single straight dislocation of length L 
and with Burger's vector b in the FLL consists of the 
core energy and of the logarithmically diverging contri- 
bution from the long-range elastic strains Eq = 
LE£)(cD + ^Tn. (L±/a)), where L± is the lateral system size, 
Eu = Kb'^ /An, K = Y^clicge is the isotropized elastic 
constant in the rescaled coordinate z — ziy^c^^/cgg, and 
CD « 1 is found numerically (044 and cge are the tilt and 
shear moduli of the vortex lattice, respectively). Bending 
of the dislocation line costs an elastic energy associated 
with its stiffness eo- Hence, the single directed disloca- 
tion line ~ parameterized by its displacement field uu(z) 
- is described by the Hamiltonian 

HdIud] =^0 + J dz^eoid^uo)^ (1) 

where the stiffness en ~ Eulnil/kza) has a logarithmic 
dispersion due to the long-range strain field. 

Let us now consider the driving force on an edge dislo- 
cation with b\ \x if a transport current j\\y is sent through 
the sample. Due to Maxwell's equation V x i? = the 
driving current creates a magnetization gradient which, 
in turn, gives rise to shear strains in the vortex lattice: 
dxUy = dxtt = a^-g- Resulting shear stresses give rise to 
a glide-component of the driving Peach-Kohler force 

F^"-^^ayxb^bKa—^ (2) 
c r> 

per length (compression stress can only give rise to dis- 
location climb which can be neglected as slow process 
requiring diffusion of interstitials ]20||). 

The displacements induced by the magnetization gra- 
dient can only be accommodated by the creation of a 
stationary superstructure of regularly spaced dislocation 
bands with Burger's vectors with an y-component ||2l[] in 
which our single "test" -dislocation is moving. The su- 
perstructure is similar to the grain boundaries formed 
in an bent atomic crystal Because such dislocation 
bands are essentially free of shear stresses [|l9| our test- 
dislocation experiences the purely elastic driving force 
everywhere in between the bands such that the presence 
of the superstructure does not affect the glide motion. 




FIG. 1. Dislocation bands obtained from numerical inte- 
gration of the time-dependent Ginzburg-Landau equation 
with current j\\y. 



In the presence of a random pinning potential Vpin (r) in 
the vortex array also the dislocation experiences Peach- 
Kohler- type pinning forces. To calculate these we have 
to determine the random stress exerted by the pinning 
potential in a frozen-in elastic displacement configura- 
tion Uei(R,z) of the vortex lattice: T^m(R + Ue/,2:) — 
afj^CR + Uei, z)'ViUei,j- The distribution of the pinning 
stresses is given by the quenched distribution of the elas- 
tic displacements Uei- This distribution is determined 
by the different spatial scaling regimes of the dislocation- 
free collectively pinned vortex array: (i) Small distances 
where vortex displacements u are smaller than the co- 
herence length ^ and perturbation theory applies [ p2[ . 
(ii) Intermediate scales where C ~ ~ ^ ^ind disorder 
potentials seen by different FLs are effectively uncorre- 
lated. This regime is captured in so-called random man- 
ifold (RM) models [^|jlo|, leading to a roughness G'(r) — 
((Ue,(r)-Ue,(0))2) ~ a2(r/i?,)2C«M ^^ere (rm ~ 1/5 
for the d — 3 dimensional RM with two displacement 
components. The crossover scale to the asymptotic be- 
haviour is the positional correlation length Ra where the 
average displacement is of the order of the FL spacing: 
M « a. (iii) The asymptotic Bragg glass regime where 
the a-periodicity of the FL array becomes important for 
the coupling to the disorder and the array is effectively 
subject to a periodic pinning potential with period a [Q. 
Here the logarithmic roughness G(r) w (o/tt)^ In (er/ Ra), 
i-e., CbtG = C(log) takes over. 

For the physics of dislocations on scales > a only 
the RM and BrG regimes are relevant. We obtain ap- 
proximately Gaussian distributed quenched stresses with 
^ = and afp(k)a!p(X) = EP»(fc)(27r)35(k + k') 
with SP"'(fc) = K'^k^G{k), i.e., 

determined by the elastic correlations G{k) with a numer- 
ical constant Brm- The RM-result holds for kRa > 1 the 
BrG-result for kRa < 1. 

To derive the correct Peach-Kohler pinning force on 
the dislocation it is crucial not only to consider the "di- 
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rect" quenched pinning stresses <Tij (r) but also the elastic 
stresses crfj themselves which are responding to the same 
pinning potential and hence tend to relax (longitudinal) 
components of the stress. A lengthy calculation shows 
that if both components are properly added the pinning 
Peach-Kohler force on a dislocation element dR 



(4) 



is rotation-free (V x dFP™)^ — such that its potential 
- the pinning Hamiltonian - can be defined: 



stability of the BrG regime with respect to dislocation 
formation The detailed stability analysis for both 

regimes was given in Ref. the BrG-AVG transition 

is weakly first-order, and dislocations proliferate with the 
density Ud — defined by the scale Ra of crossover be- 
tween the unstable RM and the stable BrG regimes. 

Now we extend our scaling analysis to the dynamic be- 
haviour of the driven dislocation line in quenched disor- 
der. Minima of the pinning potential are typically dis- 
tances ud ^ CL, given by the core size of the disloca- 
tion, apart. A dislocation confined within one minimum 
MD < a is collectively pinned on a plastic pinning length 



airWkH'i-^) - k-'YF'^{k)5uk'5iv = K'G{k)5kk'5iv (5) 

The free energy Hd [vld] [u^,] - / dzY'^""^ ■ vld from 

(0); ^'^d (H) finally provides an adequate description 
for the study of an ensemble of pinned dislocation lines. 

We start with statics and discuss roughening of the 
dislocation line in the presence of disorder. The typical 
pinning energy fluctuations upon displacing a dislocation 
segment L over a distance ud (in the xz glide plane) can 
be calculated from 



L 



pi 



Ra 

a 



<Ra 



(9) 



L U£ 







El^^{L,UD)^b^K'LuD Jdz Jdx yA^fc2g(fc)e*'=^-+*'=^- 

(6) 



2. JRM: iun/Raf^''" 



whereas the corresponding elastic bending energy of the 
dislocation is Eei{L, ujj) ~ ^Ejj In {L / a)uj~, / L (|^). Opti- 
mization gives a dislocation roughness 



RM: {L/Ra 
BrG: In 



-2/3 



|2Cf;M /(3 — 2Cjim) 

(L/a) 



O(log) 



(7) 



i.e., exponents (d ~ for RM scaling {L < Ra) 
and Cd ~ 1 ~ log'^/^ for BrG scaling. An instabil- 
ity with respect to dislocation proliferation is signaled 
by anomalous energy gains if C-D > 1, i-e., in the RM- 
regime. In the BrG-regime the energy balance is more 
subtle and to conclude C-D < 1 or stability one can 
convert the result (|^) into an approximate renormaliza- 
tion (RG) scheme: the energy gain due to roughening is 
A£; - EdLItT^/^ [L/a) - EoLioiL)-^/^ , where the 
logarithmic correction is identical to the dimensionless 
line tension cd = £d/Ed on the scale L. Interpret- 
ing AE / {EdL) as disorder correction to the line tension 
ioiL) on the scale L and summing these corrections suc- 
cessively on each scale, together with the bare tension 
e%{L) = In {L/a), one obtains an integral RG equation 



inilnL) = f d^(l±ei5(^)-i/3) 



(8) 



equivalent to the result of Ref. ||T^. Integration shows 
that corrections to i%{L) are irrelevant and hence the 



which is the analogue of the Larkin pinning length Lc — 
of the single vortex where S is the dimension- 
less pinning strength of Ref. To depin the disloca- 
tion the driving force (j^) has to exceed the pinning force 
FP™(Lp/) ~ ^Lpi {a/Ra)^'^""'^^ on a segment of length 
Lpi . This determines a critical plastic current jpi 



Jpi 



c B 
Btt^ a \Ra 



4Chm/3 



Jo 



f2TTH, 



c2 



-7/4 



<55/i« (10) 



where jo — gJi^^ -ffc2/C is the depairing current. Com- 
paring this result to the depinning current jc — jo^^^^ 
for a single vortex, one finds jpi < jc for typical disor- 
der strengths S ~ 10^'^. This shows that plastic motion 
of dislocations will set in even if the vortices themselves 
are still pinned. Thus plastic motion is the dominating 
transport mechanism in the dislocated AVG phase. 

However, the existence of a current jpi > also demon- 
strates that dislocations can be pinned at low currents and 
plastic motion for j <C jpi will only occur through activa- 
tion over diverging plastic energy barriers Upi{j) ^ j^^''^ 
giving rise to plastic creep. The typical segment size L{j) 
for activated motion at j <C jpi is determined by bal- 
ancing the energy gain by the driving force pdnvej^i+Co 
against the pinning energy Epin{L,UD{L)) ~ EdL'^'^"~^ 
of the rough dislocation line. This yields 



Upi{]) ^EDa[ — 



2Chm/3 



f) 



(2Cd-1)/(2-Ci5) 



(11) 



and we obtain the scafing law /Xp/ = (2(^_d — l)/(2 — ^d) (as 
for single vortex creep Q) relating the plastic creep expo- 
nent fipi to the dislocation roughness. We find /ipi = 
in the RM-regime {L{j) < Ra) and fipi = 1 in the BrG- 
regime {L{j) > Ra). Both exponents are considerably 
larger than their counterpart ^ = 2/11 for elastic single 
vortex creep showing that plastic creep rates are much 
smaller than elastic creep rates. 

So far we focused on a single dislocation. Now we 
turn to interacting dislocations. On large scales exceeding 
the dislocation spacing Ro, which varies from Rjj ~ Ra 
at the AVG-BrG transition to Rd ~ a at the critical 
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point H^] , interactions become essential and plastic creep 
is governed by motion of dislocation bundles in a glide 
plane (xz-plane). Deformations ud{x,z) of such a 2D 
bundle can be described by an elastic Hamiltonian with 
tilt modulus ~ Ejy/Rjy and the compression modulus 
Kx ~ R£)d'^^{Rjjf{RD)) which can be calculated from 
the dislocation free energy /(Rd), see Ref. Q (in the 
absence of disorder one finds ~ Ejj/Rd). Including 
the pinning energies we obtain the Hamiltonian 

n[uD{x, z)] = I dxdz^ {K^id^uo)^ + K,{d,UD)^} 



E 



dznp;''[h„UD{iRD,z)] (12) 



The dislocation bundle contains dislocations of opposite 
signs with the same density to avoid the accumulation 
of stress. Hence the sum over the dislocation index i 
in ( p^ ) goes over alternating Burger's vectors bi||x. On 
scales Lx ^ Rd dislocations couple effectively as dipoles 
to disorder and wc obtain for the bundle disorder en- 
ergy fluctuations E^^^{Lz, LxtUd) ^ Epm{Lz,UD) with 
Epin{Lz,UD) from (^). This has to be balanced against 
the elastic energy Eei{Lz, Lx,ud) — ^/ K^KyU^ with 
Lz ~ Kzl K^Lx resulting in a roughness 



2C^ 



un{Lz)-L\/^R^' 



RM: 
BrG: 1 



1/3^2/3 \ 3-2CRM 



(13) 
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for 



i.e., a reduced dislocation bundle roughness Qd 
RM scahng (L < Rg) and C-D ~ ^ for BrG scahng. It 
is easy to establish the corresponding plastic creep ex- 
ponents for bundle creep along the same lines as for the 
single dislocation, in particular one finds the same scaling 
relation = {2(^o d — 2) / [2 — C,d) as for d-dimensional 
vortex bundles §. This gives ^pi = ^ in the RM- 
regime and /ip/ = | in the BrG-regime. The crossover 
from single dislocation to bundle scaling happens for cur- 
rents j < where L{j^) ~ ^jK^jK-j-Ru- For Rd ^ a 
one finds ~ jpi, which means that only plastic bundle 
creep can be measured above the critical point defined 
by Rd — OL- Though we do think that our results apply 
rather to the YBCO-compound than the strongly layered 
BSCCO-compound, in Ref. a creep exponent strik- 
ingly close to fipi = I has been measured within the AVG 
phase. 

In the high-field regime above the critical point there 
is no melting phase transition but a crossover in the re- 
sistive behavior due to thermal depinning of dislocations 
will still be observable. The single pinned dislocation line 
thermally depins if {u]j)th{Lpi) ~ a?, i.e., when the ther- [17: 
mal fluctuations of a segment of the size of the plastic 
pinning length increase beyond a lattice spacing which 
leads to considerable smoothing of the pinning poten- 
tial. This happens at a plastic depinning temperature 
Tpi ~ EDa{a/Raf^'^. Using results of Ref. Q, we 
find Tpi ~ Tm — I.&Edo, at the critical point, this sug- 
gests that at the critical point the first-order melting line 



terminates and transforms into a thermal depinning line 
Tpi{B) ^ B~^/^ describing the dynamic crossover to re- 
duced depinning currents and thermal activation over all 
plastic barriers Upi{j) < T. 

In conclusion, we have developed a theory of plastic 
creep in terms of the actual dislocation dynamics in the 
pinned vortex lattice. We find infinite barriers for plastic 
vortex transport in agreement with experimentally ob- 
served low creep rates or apparently high critical cur- 
rents. The obtained results are relevant for and are easily 
extended to other systems where glassy dynamics is con- 
trolled by topological defects, for example, charge density 
waves in disordered crystals and/or work-hardened solids. 
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